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Abstract 


This paper is the first part of a two-part paper and presents the development, calibration and validation of a two-dimensional isothermal mechanistic 
model of a composite yttria/scandia-stabilized zirconia anode-supported multiple layers solid oxide fuel cell (Ni-YSZ|Ni-ScSZ|ScSZ|LSM-ScSZ). 
This model was developed to describe the intricate interdependency among the ionic conduction, electronic conduction, multi-component species 
transport, electrochemical reaction processes and electrode microstructure for intermediate temperatures operation between 750 and 850°C. This 
model takes into account the fact that the electrochemical reactions take place throughout the electrodes and not only at the electrolyte/electrode 
boundaries. The model was calibrated using experimental polarization curves and then validated by comparing each cell component polarizations 
(anodic, cathodic and electrolyte) determined from the simulation and from specific experiments using a symmetric cell and EIS measurements. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


A solid oxide fuel cell (SOFC) is an energy conversion device 
that converts the chemical energy of a fuel directly into elec- 
tricity [1]. With rising fuel prices and stricter emission control 
regulations, solid oxide fuel cells become even more attractive 
due to their high efficiency, low environmental impacts and fuel 
flexibility. 

Electrochemical reactions in SOFC are taking place in a 
PEN (positive electrode|electrolyte|negative electrode). In order 
to overcome the problems associated with high temperature 
operation of state-of-the-art SOFC (around 1273 K), a growing 
number of researchers are focusing their efforts on intermedi- 
ate temperature solid oxide fuel cells (IT-SOFCs) for operation 
between 823 and 1073 K. IT-SOFCs allow for a wider range 
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of materials and more cost-effective SOFC fabrication meth- 
ods [2,3], However, the ionic conductivity of the electrolyte 
decreases as the operating temperature is lowered. The higher 
ohmic overpotential at lower operating temperatures could be 
reduced by using composite ceramic electrolytes that have 
higher ionic conductivity than conventional YSZ electrolyte, or 
by adopting an electrode-supported configuration with a thinner 
electrolyte. One promising configuration for IT-SOFCs is an 
anode-supported SOFC where a thin electrolyte of thickness in 
the ranges of 8—15 um is deposited on a thick anode [4]. In addi- 
tion, the ionic conductivity and anodic reaction activity could 
be further improved by adopting a scandia-stabilized zirconia 
(ScSZ) electrolyte and a Ni/ScSZ composite anode [5,6]. The 
operation of SOFC involves complex chemical, electrochemical 
and mass transport processes. The cell performance is strongly 
affected by several irreversible losses including ohmic losses 
due to ionic and electronic charge transfer resistances, activa- 
tion losses due to irreversibility of electrochemical reactions at 
the three-phase boundary (TPB), and concentration losses due 
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Nomenclature 

c concentration (mol m~?) 

D diffusion coefficient (m? s7!) 

E activation energy (J mol!) 

E? standard cell potential (V) 

F Faraday constant (96,485 C mol~!) 

i current density (A m~?) 

itrans transfer current density (A m~?) 

io exchange current density (A m~?) 

I current (A) 

Mi molecular weight of species i 

Ne number of electrons participating in the reaction 

N molar flux (mol m~? s7!) 

p pressure (Pa) 

Q source term for charge balance equations (A m~*) 

r average pore size (wm) 

R cell radius (m) 

R gas constant (8.314 J mol`! K7!) 

Ri source term for mass balance equations 
(kg m~? s7!) 

S cell area (m~?) 

STPB TPB active area per unit volume (m? m~?) 

T temperature (K) 

V voltage (V) 

Wi mass fraction of species i 

Xi mole fraction of species i 


Greek letters 


a transfer coefficient 

B adjustable parameter in exchange current density 
formulations (see Eqs. (1) and (2)) (Q27! m~?) 

E porosity 

n overpotential (V) 

p density (kg m~?) 

o conductivity (S m7!) 

T tortuosity 

Subscripts 

ac anode chamber 

act active layer 

an anode 

avg average 

ca cathode 

cc cathode chamber 

elec electronic 

electrolyte electrolyte 

inter interface 

ion ionic 

leak leak 

op open circuit 

p pore 

ref reference 

sp support layer 


symmetry axis of symmetry 


theoretical 


to mass transport resistance in the electrode, especially for thick 
anodes as in an anode-supported SOFC [7]. In reality, leak over- 
potential associated with parasitic loss due to current leakage, 
gas crossover, and unwanted side reactions must also be consid- 
ered. High SOFC performance relies on optimal electrochemical 
reactions and mass transport processes. Since experimental stud- 
ies on SOFC are expensive, time-consuming and labor-intensive, 
quantitative mechanistic models for the cell PEN structure are 
essential for SOFC technology development. A validated mech- 
anistic model offers not only a means to better understand the 
complex physical phenomena governing fuel cell performance 
that are not readily accessible experimentally; they also can be 
useful design tools. 

In recent years, several SOFC models have been developed to 
study reactions and transport phenomena within the PEN struc- 
ture. These models differ widely in terms of their complexity 
and comprehensiveness. 

Tanner et al. [8] developed an empirical model based on the 
works of McDougall [9] to calculate the effective charge transfer 
resistance of an electrode as a function of intrinsic charge trans- 
fer resistance, ionic conductivity of the electrolyte and electrode 
thickness. In order to predict the concentration polarization in 
detail, Virkar et al. [10] and Kim et al. [11] considered porous 
media gas phase transport processes based on Fick’s model to 
predict concentration overpotentials, and used Tanner’s method 
to calculate activation overpotentials. Although those empirical 
models could predict the cell performance conveniently they are 
limited to particular cells and cannot describe in detail the effects 
of operating and design parameters. 

One-dimensional PEN models adopt conservation laws to 
describe the reaction and transport processes occurring within 
the cell and could be used to identify the effects of each layer 
in the PEN. Costamagna [12,13] developed a one-dimensional 
continuum micro-scale models in which the effective proper- 
ties were estimated based on statistical properties of random 
packing systems of binary spherical particles and percolation 
theory. Chan et al. [14,15] adapted this model to consider the 
effects of concentration overpotential by taking the molecular 
and Knudsen diffusion into account. These models were used 
for detailed simulation of only one electrode (anode or cath- 
ode). These continuum micro-scale models were extended to 
the whole PEN by Nam and Jeon [16] and Shi and Cai [17] to 
build a bridge connecting micro-scale to macro-scale models. 
However, a one-dimensional PEN model could only represent 
the parameter variation along the direction of the PEN thickness. 
For instance, they could not be used to calculate the distribution 
of the electric current path inside the PEN. In addition, such 
one-dimensional models cannot model the effects of asymmetric 
electrodes, which were thus omitted. Therefore, a comprehen- 
sive multi-dimensional model is necessary when attempting to 
simulate the intricate interdependency among the ionic con- 
duction, electronic conduction, gas transport phenomena, and 
electrochemical processes. Numerous multi-dimensional SOFC 
models exist in the open literature with varying sets of assump- 
tions depending on the objectives of the simulation [18-23]. 
However, in order to simplify the calculation, these published 
multi-dimensional models treated the reaction zone as boundary 
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conditions. This is not consistent with the fact that the reaction 
zone layer in an SOFC is spread into the electrode to some 
distance from the electrode/electrolyte interface [14,15,17,24]. 
In addition, since the electrode microstructure has a direct 
impact on the cell performance [12], including micro struc- 
tural characteristics in the model can provide insights that relate 
microstructure to performance. Such a model would be partic- 
ularly useful for cell design and optimization. Unfortunately, 
few of the existing multi-dimensional PEN models relate the 
electrode microstructure to the cell electrochemical reaction 
characteristics. 

Furthermore, for a given mathematic model, although most 
of the model parameters are determined from experiments or 
from literature, others are only estimates. Some of the parame- 
ters are adjusted to ensure good agreement between the model 
results and the experimental data, for examples exchange cur- 
rent density and tortuosity [17]. The uncertainty in the estimated 
parameter will surely influence the model reliability and accu- 
racy. The more the parameters are determined independently 
through experiments or characterization techniques, the more 
robust the model. 

In this paper, an isothermal mechanistic model that considers 
the cell microstructure and reactions to occur over the whole 
electrode was calibrated and validated using experimental 
data obtained on a button cell consisting of Ni-YSZ|Ni- 
ScSZ|ScSZ|LSM-ScSZ multiple layers. This model has 
already been presented in detail in a previous paper [25]. 
However, contrary to the model presented earlier on, the present 
model can accommodate multi-components. In addition, an 
important feature of the validation process presented here is that 
the model is not only validated using the overall polarization 
curve as it is usually done in the literature, but it is also validated 
by considering overall losses over each cell component, namely 
anode, cathode and electrolyte, and thereafter referred to as 
anodic, cathodic and electrolyte overpotentials, respectively. 
Note that in the “electrolyte” overpotential we also include 
the loss due to contact resistance between the electrodes and 
electrolyte. The methods to isolate each cell component over- 


O, 


TPB 


Cathode layer (LSM/ScSZ) 
Electrolyte layer (Dense ScSZ ) 


Anode active interlayer (Ni/ScSZ) 
Anode support layer (Ni/Y SZ) 


potential from experiments and from simulation are described 
in Sections 5.2 and 5.3, respectively. 


2. PEN structure and experimental apparatus 
2.1. Anode-supported SOFC PEN structure 


Fig. | illustrates the PEN structure of the anode-supported 
button cell. The PEN used in this paper consists of a Ni/YSZ 
anode support layer (680 um), a Ni/ScSZ anode active inter- 
layer (15 wm), a ScSZ thin-film electrolyte layer (20 wm) and a 
LSM/ScSZ cathode layer (15 um). Except for the cathode layer 
which has a diameter of 1.4 cm, all other layers have a diameter 
of 2.6cm. The electrodes are composed of electronic and ionic 
conductor particles and are porous to facilitate the transport of 
fuel and oxidant from the gas channel to the three phase bound- 
aries where the electrochemical reactions occur. The electrolyte 
is dense to keep the air and fuel gases separated and thus allow- 
ing for an oxygen concentration difference between the anode 
and the cathode. Oxygen ions are produced at the TPB sites 
near the cathode/electrolyte interface, and are transported by a 
solid-state migration mechanism through the electrolyte to the 
anode/electrolyte interface, where oxygen ions react with the 
fuel. Product molecules are then transported back to the fuel 
channel through the pores. 


2.2. Cell preparation 


Commercial NiO (Inco Ltd., Canada) was used with 8 mol% 
yttria-stabilized ZrO2 (YSZ, TOSOH, Japan) powder for prepar- 
ing the support layer and with scandia-stabilized zirconia, 
Zro.89Sc0.1 Ce0.0102—x (ScSZ, Daiichi Kigenso Kagaku Kogyo, 
Japan), for the anode active layer. The ratio of the mixtures 
(NiO-YSZ and NiO-ScSZ) was 50 wt% NiO and 50 wt%- 
stabilized zirconia. The ScSZ powder was also used to prepare 
the electrolyte film. 

The slurries for the tape casting process were prepared by a 
ball milling process that included two steps. In the first step, all 


Three phase 
boundary 


See 33] electronic conductor 


ionic conductor 


Fig. 1. Schematic representation of the anode-supported SOFC button cell used in this study. 
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the above-mentioned ceramic powders were homogenized in a 
planetary mill for 2 h with dispersant in a mixture of methy] ethyl 
ketone (MEK) and ethanol (EtOH). To form sufficient porosity 
in the anode support and active layers, a certain amount of rice 
starch was added as pore former to the mixture of NiO-Y SZ and 
NiO-ScSZ. Other organic additives, such as polyvinyl butyral 
(PVB) and a mixture of polyethylene glycol (PEG) and dibutyl 
o-phthalate (DOP) were added as binder and plasticizer, respec- 
tively, in adequate quantity and proportion, and then milled for 
another 2h. 

Prior to tape casting, the slurries were vacuum pumped for 
about 2 min in order to remove any remaining air. The ScSZ film 
was casted first onto the glass plate by a “doctor blade” method 
and allowed to dry in air for several minutes; the anode functional 
and active layers were prepared similarly. After drying overnight 
at room temperature, the multilayer green tape was detached and 
co-sintered at 1400 °C in air for 4h. 

The LSM (Inframat advanced materials, USA) and ScSZ 
powders in a mass ratio of 50:50 were mixed in a planetary mill 
with ethanol for 3 h to ensure random distribution of each phase. 
The dried mixture was subsequently grounded in an agate mor- 
tar with ethyl cellulose and terpineol to prepare the paste, which 
was screen-printed onto the sintered ScSZ layer and sintered at 
1200°C for 3h to form the cathode. 


2.3. Test setup 


Fig. 2 shows the schematic representation of the test station 
used for evaluating the performance of the button cell test. 


Alumina tube —— 


The cell was located between two alumina tubes. Pt wires 
were used as voltage and current probe. The Pt mesh was used 
as the cathode current collector and was fixed to the porous 
cathode with platinum paste. The porous Ni felt was used as 
the anode current collector and was fixed to the anode support 
layer with platinum paste. Before testing, the anode was fully 
reduced at 850°C in a 20:80 H2:N2 mixture (total flowrate of 
50sccm). During the actual tests, H2 mixtures were humidi- 
fied at room temperature (30°C, 4% H20) and used as fuel; 
oxygen or air was used as oxidant. The flow rates of fuel and 
oxidant were both kept at 50 sccm. The polarization curves were 
measured using four-probe method with an Electrochemical 
Workstation (IM6e, Zahner-Elektrik GmbH, Germany) for tem- 
peratures ranging from 700 to 850°C. The measurements were 
initiated after the system was stable under a constant voltage 
of 0.7 V for 30 min. Electrochemical impedance spectroscopy 
(EIS) was performed using an amplitude of 10 mV over the fre- 
quency range 0.1 Hz to 100 kHz. The ohmic resistance of the 
whole cell was estimated from the high frequency intercept of 
the impedance curves. 

After the cell tests, the cell surface and cross-section, as well 
as the elemental distribution, were characterized by scanning 
electronic microscope (SEM) and energy dispersive spectrom- 
eter (EDS) using an electron probe microanalyzer (JSM-6460, 
JEOL, Japan). The actual thickness of the cell components were 
also determined using SEM. The porosity and tortuosity of the 
different layers before and after the cell tests were determined 
using mercury porosimetry (AutoPore IV 9510, Micromeritics, 
USA) and image processing. 
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Fig. 2. Schematic representation of the test setup. 
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3. Model development 
3.1. Model assumptions and model geometry 
The assumptions made in this model are the following: 


(1) The PEN operates at steady state. 

(2) All gas mixtures are considered as ideal gases. 

(3) The reaction active sites are assumed to be uniformly dis- 
tributed in each electrode layer. The two conducting phases 
(electronic and ionic) are considered to be continuous and 
homogeneous in each layer. 

(4) The temperature is uniform in the PEN, thus the model is 
assumed isothermal and all physical properties are evaluated 
at a fixed cell temperature. 

(5) Convection flux is neglected in the porous electrode com- 
pared to diffusion. Pressure gradients in the porous electrode 
are also neglected. 

(6) The boundary conditions for potential and species concen- 
trations are assumed uniform at the electrode/gas channel 
interface. 


In addition, it should be noted that the reaction zone may 
not be restricted to the anode active interlayer but may spread 
into the anode support layer, especially if the active interlayer is 
relatively thin. Thus, the governing equations for the anode- 
supported layer are the same as for the anode active layer, 
although some parameters may be different. 

Due to symmetry, a two-dimensional axial symmetry coor- 
dinate is adopted, as shown in Fig. 3. In this figure the labeling 
of the calculation domains and boundaries are also shown. The 
type of boundary conditions for each domain will be described 
later in the paper. 

With the above assumptions and model geometry, the SOFC 
PEN model is formulated using charge balances, mass balances 
together with electrochemical reaction kinetics, as described in 
the following sections. 


3.2. Governing equations 


Most of the governing equations have been presented in 
detail in a previous paper [25]. Only the differences between 


the present model and the model used in [25], as well as some 
additional information are highlighted here: 


e Inthe present model surface diffusion has not been taken into 
account. 

e Material balances in the present model accommodate more 
than two components using the dusty gas model [28]. 

e The expression of anodic exchange current density (both 
anode support layer and anode active layer) is kept the same 
as in Ref. [25]: 


. Ban RT / cm, Eact,an 0.133 
i = exp |— ? an)” 1 
0,an 3F a p RT (Poo, an) (1) 


where ig,an is the anodic exchange current density, Cref,H, 
is the hydrogen concentration in the reference case and 
PO>,an is the oxygen partial pressure at the anode. For the 
support layer we kept the values found in Ref. [27], that 
is Eact,an = 120,000 J mol“! and Ban =6.17 x 10!! Q7! m-?. 
For the active layer, Eact,an is kept constant at 120,000 J mol! 
and Ban is treated as an adjustable parameter to fit the exper- 
imental data. 

e The cathodic exchange current density is expressed in a dif- 
ferent form than in Ref. [25]. It is determined using the 
expression found in Costamagna et al. [27]: 


BeaRT Enact, 
aF exp acboa ) Cpo ca)?” (2) 


10,ca = 


where io,ca is the cathodic exchange current density and po, an 
is the oxygen partial pressure at the cathode. Eactca is kept 
constant at 130,000 J mol—! while Bea is also treated as an 
adjustable parameter, similar to Ban. 

e 7 is the local overpotential defined as 


n = Velec — Vion — Vref (3) 


where Vier is the relative potential difference between the 
electronic and ionic conductors at a reference state. 

e Open-circuit state was chosen as the reference state here. By 
setting Vref,an to zero, the cathode reference potential Vref.ca 
will be the open circuit voltage. However, it should be noted 
that the experimental open circuit voltage values are typically 
slightly less than the theoretical values. This difference is due 
to slight gas leaks in a single cell test and thus, this potential 
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Fig. 3. Two-dimensional axisymetrical model geometry. 
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difference is called “leak overpotential’”, denoted as njeak and 
kept constant in this study. The cathode reference potential 
can thus be formulated as 


Vref,ca = Vop + Vref,an = Vop,th — Mleak 
URA 
=0 
0.5 
— E? RT In PO3,caPH2.an 
2F PH20,an 


Neak (4) 


where Vop denotes the actual cell open circuit voltage, Vop,tn 
is the theoretical open circuit voltage, E? is the standard cell 
potential, and py, an and PH,O,an are the partial pressures of 
hydrogen and water at the anode, respectively. 


The governing equations for the electrode charge balances 
are summarized in the following equations (see nomenclature 
for the definition of each term and Ref. [25] for more details). 


e Ionic charge at the cathode: 


—V (6% ca V Vion,ca) 
= Qion,ca = —i0,ca STPB,ca 
TPB 
CO, ane F(Velec,ca = Vion,ca = Vref,ca) 
x bulk ex RT 
cO 
a -— a)Ne F(Velec,ca = Vion,ca — Vref,ca) 
exp 
RT 
(5) 
e Electronic charge at the cathode: 
-V (of, caV Velec,ca) = Qelec,ca = —Qion,ca (6) 
e Ionic charge at the anode: 
eff A 
= V(Gion.an V Vion,an) 
= Qion,an = io,an STPB,an 
TPB 
CH, ane F(Velec,an — Vion,an — Vief,an) 
x | bulk RT 
H2 
TPB 
CHO ( (1—@)ne F(Velec,an— Vion,an = fetan ) 
bulk 
HÔ RT 
(7) 
e Electronic charge at the anode: 
-V (oan V Velec,an) = Qelec,an = — Qion,an (8) 


The governing equations for the electrode mass balances 

are summarized as follows: 

e Mass balance at the anode: 
n b 

—ltrans,an STPB,an My, 


—V pwn, X Dn, jVxj = ZF (9) 


j=l 


Table 1 
Boundary conditions 


Boundary Ionic charge Electronic charge Mass balance 
IQ calec Insulation Veell,ca WO, bulk,ca> 
WN), bulk,ca 
3 NRelec/ce Insulation Insulation Insulation 
3NR ea/elec Continuity Insulation Insulation 
OD declan act Continuity Insulation Insulation 
O92. an_act/an_sp Continuity Continuity Continuity 
3Nan-sp/fe Insulation Veell,an WH3,bulk,an> 
WH20,bulk,an> 
WN3,bulk,an 
IN symmetry Symmetry Symmetry Symmetry 
Others Insulation Insulation Insulation 


n 


= —itrans,an STPB,an MH30 
—V PWHO y Dp,o, jV xj = 


= 2F 
(10) 
e Mass balance at the cathode: 
n # 
z — sanS Mo, 
-y pwo, X Do, ;Vxj _ ee Oo (11) 


3.3. Boundary conditions 


In order to solve the system of coupled partial differential 
equations for charge and mass balances, the boundary conditions 
of all outer interfaces are specified in Table | (refer also to Fig. 3). 

The difference between Veelica and Veellan is the cell volt- 
age used in the calculation. Here, we chose Veetlan =0. Wi bulk 
in the table is the mass fraction of species i in the fuel/air 
channel. 

The boundary conditions “Insulation” and “Symmetry” both 
mean that the partial derivative of the variable at the boundary 
is zero. The difference is that for “Symmetry” there is a flux 
through the boundary, while for “Insulation” there is none. 


3.4. Solution method 


The problem was solved for a given cell voltage. The outputs 
of the model are the distributions of current density and species 
concentrations. The calculations were performed using the finite 
element commercial software COMSOL MULTIPHSICS®, 
Version 3.2. 

From the current density distribution, the average current 
density was calculated as 


ing = “= = f aai (12) 
where ijocaj is the local current density and R is the radius of the 
interface based here on the cathode area. 

In order to generate a complete polarization curve, the calcu- 
lation was performed over a range of cell voltages, for which a 
corresponding average current density was determined. 
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4. Model calibration 


The process of model calibration involves the following two 
steps: (1) evaluation of model parameters such as microstruc- 
ture properties and kinetic parameters; (2) model calibration for 
a “base case” (here humidified H2 with air at 1 atm for temper- 
atures between 750 and 850°C). Ideally, all model parameters 
should be determined independently from literature data or from 
characterization techniques. Unfortunately, this is not possible 
for some parameters (e.g. leak overpotential, an intrinsic param- 
eter of the experimental setup), which then must be estimated by 
tuning their value to fit the experimental data during the calibra- 
tion step. However, once all model parameters are determined, 
they were not changed when simulating other cases as in the 
final step of the model validation procedure. 


4.1. Model parameters 


Except for the exchange current density, which expression 
was found in Costamagna et al. [26], all other kinetic param- 
eters were taken from the work of Nagata et al. [27]. This 
section describes then how the model parameters related to 
the microstructure of the different cell layers were estimated. 
Those parameters are: layer thickness, ionic conductivities and, 
in the case of the electrodes, pore size, tortuosity and electronic 
conductivities. 

A picture showing the appearance and dimension of the but- 
ton cell used in this study is shown in Fig. 4. The thickness of 
each layer was determined from SEM micrographs of the cell 
cross-section, as shown in Fig. 5. 

The pore size and porosity of the cell were characterized using 
mercury porosimetry. The mean pore diameter and porosity were 
found to be 0.46 um and 0.364, respectively. It is very difficult to 
differentiate the pore size and porosity of anode active layer and 
cathode layer from that of the anode support layer since they 
are very thin compared to the anode support layer. Actually, 
the values of pore size and porosity determined from mercury 
porosimetry are essentially characteristics of the anode support 


Fig. 4. Photograph of the button cell used in the experiment. Left: cathode side, 
right: anode side. 


(a) 


(b) 


Fig. 5. SEM micrographs after tests: (a) overall cell cross-section and (b) 
detailed view of cell cross-section. 


layer. Thus, an image processing software (ImageJ V1.34) was 
used to determine the pore size and porosity of each layer from 
SEM images based on quantitative stereology [28,29]. It should 
be noted that the porosity and pore size from image analysis 
were only used to determine their relative value compared to 
that of the anode support layer. The values used in the model 
were adjusted accordingly, based on the values of pore size and 
porosity of the anode support layer determined from mercury 
porosimetry. 

At first, the SEM images were contrasted and filtered by the 
ImageJ software. After thresholding, the binarized images were 
used to determine the porosity and the pore size distributions. 
According to stereological assumptions, the calculated porosity 
is equal to the pore volume proportion [28]. An “equivalent pore 
radius” [29] was used to quantify the pore size since the shape 
of the pores varies significantly. To reduce measuring errors, 
measurements were performed five times at different position of 
each layer. Fig. 6 illustrates a comparison between pore size dis- 
tribution obtained by image analysis and mercury porosimetry. 
The porosity, average pore area and equivalent radius of anode 
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Table 2 
Porosity, average pore area and equivalent radius evaluated by image analysis 


Sample Anode support layer Anode active layer Cathode layer 
number 
Porosity Average pore Equivalent pore Porosity Average pore Equivalent pore Porosity | Average pore Equivalent pore 
area (m2) radius (um) area (m?) radius (um) area (m?) radius (um) 

1 0.265 6.92E—13 0.47 0.261 2.83E-13 0.30 0.282 4.94E—13 0.40 

2 0.291 8.31E—13 0.51 0.277 2.71E—13 0.29 0.274 4.25E—13 0.37 

3 0.248 6.38E—13 0.45 0.271 3.27E— 13 0.32 0.272 5.79E— 13 0.43 

4 0.299 8.29E—13 0.51 0.278 2.77E—13 0.30 0.252 5.27E—13 0.41 

5 0.267 6.58E—13 0.46 0.291 3.39E— 13 0.33 0.253 5.39E-13 0.41 

Average 0.274 7.30E—13 0.48 0.276 2.99E— 13 0.31 0.267 5.13E-—13 0.40 


active layer, anode support layer and cathode layer evaluated 
from image processing are shown in Table 2. 

The results show that the pore size distribution of the anode 
support layer obtained from image analysis is somewhat differ- 
ent from that obtained from mercury porosimetry. The porosity 
determined from image analysis (ca. 0.27) is lower than that from 
mercury porosimetry (0.365). This is not too much an issue here, 
because we are only interested in knowing the relative values of 
porosity and pore size of the anode active layer and cathode 
layer, compared to that of the anode support layer. 

According to Table 2, the porosity of the three layers can be 
assumed identical. The average pore radius of the cathode layer 
(0.40 um) is considered close enough to the value of that of the 
anode support layer (0.48 um) that we considered them to be 
the same. However, the average pore radius of the anode active 
layer (0.31 um) is assumed to be 1.5 times smaller than that of 
the anode support layer. 

Particles in the electrode consist of agglomerates of differ- 
ent sizes before sintering. To simplify the calculation, the mean 
particle radius of the ionic conductor is assumed to be equal to 
that of the electronic conductor within the same layer [12]. The 
relationship between the mean particle radius, rejec, and mean 
pore size, rp, is as follows [16]: 


3(l— € 
i (13) 
2e 


ZZ/ Image analysis 
—=— Mercury porosimetry 


0.0 0.2 0.4 0.6 0.8 1.0 
Pore diameters/um 


Fig. 6. Pore size distribution of anode support layer obtained by image analysis 
and mercury porosimetry. 


Table 3 
Material conductivities 


Material Conductivity (S m7!) 


Tonic conductor 


ScSZ 6.92E4 exp(—9681/T)* 

YSZ 3.34E4 exp(—10300/7) [30] 
Electronic conductor 

LSM 4.2E7/T exp(—1150/T) [30] 

Ni 3.27E6 — 1065.37[16] 


Equivalent ionic conductivity 0.002T — 1.4483* 


of electrolyte layer 


è Curve fitting according to the experimental data. 


Table 3 lists the conductivities of ionic conductor, electronic 
conductor as well as the equivalent conductivity of the elec- 
trolyte layer (i.e. when the whole conductivity of the cell is 
reduced to that of the electrolyte layer), which was determined 
by EIS spectra at open circuit state. 


4.2. Model calibration 


During model calibration, a few parameters have been 
allowed to vary and be tuned to best fit the experimental data 
under one set of operating conditions referred to as “base case”. 
The operating conditions of the base case are listed in Table 4. 
The parameters to be tuned were the leak overpotential, cathode 
and anode tortuosities, and kinetics related parameters since the 
material used in the active anode layer is different from that of 
Costamagna and Honegger [26] and Nagata et al. [27], whose 
kinetic parameters were used as a starting point. The model 
parameters tuned in the base case are summarized in Table 5. 
The cell performance at three temperatures (750/800/850 °C) 
were calculated and measured experimentally while the fuel was 
humidified at a constant temperature of 30°C. 


Table 4 

Operation conditions in the base case 

Parameters Value 
Pressure, p (Pa) 101,325 
Temperature, T (°C) 750/800/850 


96% Hz and 4% H20 (saturated 
hydrogen at 30°C) 
21% Oz and 79% No 


Fuel composition 


Oxidant composition 
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Table 5 

Tuned model parameters 

Parameter Value 
Leak overpotential, nieak (V) 0.03 
Ban in Eq. (1) (Q7! m~?) 6.8E12 
Bea in Eq. (2) (Q7! m~?) 5.8E10 
Anode tortuosity 16.5 
Cathode tortuosity 3 


o 850°C,Experimental; —"— 850°C, Simulated 
4 A 800°C,Experimental; —4— 800°C, Simulated 
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Fig. 7. Comparison between modeling and experimental data for the base case. 


The polarization curves shown in Fig. 7 show that the experi- 
mental and modeling results at different temperatures agree well. 
Obviously, the operating temperature significantly affects the 
cell performance. Although the open circuit voltage decreases 
when increasing temperature, cell performance at higher temper- 
atures is better because not only are the electrochemical reactions 
and mass transport faster, but also the ionic conductivity of the 
electrolyte is considerably higher. 


5. Model validation 


The experimental polarization curves represent only the gross 
behavior of each component of the single button cell. The 
model was therefore validated by comparing each cell compo- 
nent polarization loss (anodic, cathodic and electrolyte) obtained 
experimentally and through simulation. Note that to determine 
experimentally the various polarizations, as explained in the next 
section, a different set of experiments using a different cell was 
used. 


5.1. Determination of each cell component polarization 
from experiments 


To determine each polarization loss experimentally, a sym- 
metric cell, shown in Fig. 8, was used to estimate first the 
cathodic polarization. The symmetric cathode was printed on a 
commercial ScSZ film in the same way as described previously 
for the single button cell. The total ohmic resistance was mea- 
sured by EIS tests on the symmetric cell from the high frequency 
intercept of the impedance curves. By subtracting the ohmic 
polarization from the measured polarization and then dividing 
by two, the cathodic polarization was determined. The anodic 


_7~— Symmetrical cathode 


_-_ Electrolyte layer 


Current probe Current probe 


Voltage 


Fig. 8. Symmetric cell used to determine the cathodic polarization. 


polarization was then estimated by subtracting the ohmic polar- 
ization, cathodic polarization and leak overpotential from the 
total polarization of the single cell. The symmetric cell used for 
polarization measurement should be carefully prepared to make 
sure that the cathode in the symmetric cell and the one in single 
button cell were identical. Through EIS measurements at high 
frequency, considering that the electrodes are good conductors, it 
is reasonable to assume that the resistance thus determined (high 
frequency intercept of the impedance curves) include essentially 
the contact and electrolyte ohmic resistances. The polarization 
due to both contact and electrolyte resistances was then cal- 
culated from the value of the resistance determined from EIS 
measurement. 


5.2. Determination of each cell component polarization 
from simulations 


The local overpotential was defined in Eq. (3). The cathodic 
and anodic overpotentials were calculated as 


Nan = \(Veleclan—sp/ac E Vion lelectrolyte/an—act) = Vief,an| (14) 
Nea = \(Veteclea/ac = Vion lelectrolyte/ca) = Vief,cal (15) 
The ohmic polarization was calculated as 


Nelectrolyte = (Vion lelectrolyte/ca = Vion |electrolyte/an_act) (16) 


Recall that, as defined in Section | and as seen in Eq. (16), the 
“electrolyte” overpotential (electrolyte) also includes the loss due 
to contact resistance between the electrodes and the electrolyte. 

The total cell overpotential was then calculated as 


Notal = Nan + Nca + Nelectrolyte + Mleak (17) 
5.3. Comparison experiments/simulations 
It can be seen from Fig. 9 that the calculated total overpoten- 


tial and the overpotential of each component all agree well with 
the experimental data in the base case at 800°C. 
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Fig. 9. Modeling and experimental results of cell overpotentials. 


As a mathematical trick, if the source terms of the anode and 
cathode mass balance equations (Eqs. (9)—(11)) are set to zero, 
the fuel and oxidant concentrations are uniform throughout the 
electrodes and equal to the concentrations at the electrode/gas 
chamber interfaces. In this way, the effects of concentration 
overpotential could be separated. The simulated detailed contri- 
butions of each PEN component to cell overpotential (including 
concentration overpotential) in the base case at 800 °C are shown 
in Fig. 10. 

It is seen from Fig. 10 that the cathodic activation overpoten- 
tial and the ohmic overpotential dominate for SOFC operation 
at 800°C with air and current density below 10,000 A m~?. The 
activation polarization of the anode is between 40 and 50% of 
the cathodic activation overpotential. The cathodic concentra- 
tion polarization in the base case can be neglected because the 
cathode is thin enough. The anodic concentration polarization is 
more important than the cathodic one because of the thick anode. 
Nonetheless, compared to the anodic, cathodic and ohmic losses, 
the anodic concentration overpotential remains small (less than 
5% of total overpotential) in the base case for our button cell. 
However, the anodic concentration overpotential increases with 
the current density. 

The effect of temperature on the cathodic overpotentials is 
shown in Fig. 11 for temperatures between 750 and 850°C. 
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Fig. 10. Contributions of each PEN component to cell overpotential (simulation 
results). 
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Fig. 11. Calculated and experimental cathodic overpotential at 750, 800 and 
850°C. 


Both calculated and experimental results show that the operat- 
ing temperature significantly affects the cathodic overpotential, 
which increases as the temperature decreases, as expected. In 
addition, this figure shows that the model is valid over a wide 
range of temperatures, but somewhat understimates the cathodic 
polarization at 850 °C, and somewhat overestimates the cathodic 
polarization at 750 and 800°C. 


6. Conclusions 


A two-dimensional isothermal mechanistic PEN model of 
an anode-supported solid oxide fuel cell was developed to 
describe the intricate interdependency among the ionic conduc- 
tion, electronic conduction, multi-component species transport, 
electrochemical reaction processes and electrode microstruc- 
ture. The model was calibrated and validated for an IT-SOFC 
button cell, comprised of Ni-YSZ|Ni-ScSZ|ScSZ|LSM-ScSZ 
multiple layers. Five parameters were tuned, namely leak 
overpotential, pre-exponential factors of anodic and cathodic 
electrochemical kinetic expressions (the B’s) and anodic and 
cathodic tortuosities. Tuning of these parameters led to a good 
match between the overall polarization curves obtained exper- 
imentally and from simulation for humidified hydrogen (4% 
H20) operating at temperatures between 750 and 850°C. The 
model was then validated by comparing each cell component 
polarizations (anodic, cathodic and electrolyte) determined from 
the simulation and from specific experiments using a symmetric 
cell and EIS measurements. This validation process showed that 
the model was able to describe accurately each polarization. In 
part II of this two-part paper further validation is carried out over 
a wider range of operating conditions (e.g. different fuel and oxi- 
dant compositions) and the model is used to further discuss the 
effect of cell microstructure on cell performance. 
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